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ABSTRACT 


This report considers longitudinal waves traveling in a cylindrical 
rod. After a review of the classical results of wave propagation in 
unbounded elastic solids and in elastic cylinders, the subject of thermal 
modification of elastic properties is undertaken. The resultant effects 
upon the purely elastic wave motion consists of a change in the propaga- 
tion velocity and the addition of a damping effect upon the mechanical 
energy of the elastic wave due to heat condution. The thermoelastic 
effects upon a longitudinal elastic wave are shown to be small in tenfts 
of the propagation velocity. However, more significantly, the thermo- 
elastic damping effect is large for very high frequency waves traveling 
in small diameter bars. The behavior of the thermoelastic damping coef- 
ficient is linked to a "thermoelastic bar number," developed herein. 
Finally, approximations are developed, based upon the value of the 
thermoelastic bar number, to predict the behavior of the thermoelastic 
damping coefficient as a function of frequency. 
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INTRODUCTION 

The purpose of this report is to ascertain the effects of thermo- 
elastic damping on the propagation of longitudinal waves in cylin- 
drical rods. 

As a wave propagates down the length of a rod, the variation of 
compress ional strain in the medium generates a small amount of heat 
flow which increases the entropy of the bar and dissipates a portion 
of the wave energy. This process is referred to as thermoelastic 
damping. In the past history of the development of solid mechanics, 
the magnitude of thermoelastic damping was considered extremely 
small. Therefore, emphasis was placed upon elastic wave propagation 
while thermal wave propagation was studied by thermodynamicists in 
terms of the heat transfer problem. Only recently, starting in the 
late 1950's, has the coupling of the thermal and elastic waves been 
considered. Thermoelastic wave propagation has been studied mainly 
for extended solids by H. Deresiewicz (1957) and by Chadwick and 
Sneddon (1958) and in a review article by Chadwick (1960), which 
also covers bounded solids. In particular, Chadwick extended the 
Pochammer-Chree (1886) analysis of longitudinal elastic wave propa- 
gation in infinite circular rods to the thermoelastic case. 

However, the Pochammer-Chree elastic solution and the Chadwick 
thermoelastic extension are unnecessarily complicated for present 
engineering applications and therefore present a problem to which 



this report is addressed, namely, an engineering understanding of 
thermoelastic wave propagation and damping. 

The approach to this problem taken in this report is to fina a ,« 

satisfactory engineering approximation for elastic waves and then 
extend that solution to incorporate the coupling due to the thermal 
wave. 

The literature abounds with engineering solutions of elastic 
wave propagation. Several of these methods are reviewed in 
Chapter II, and a particularly simple one-dimensional solution by 
Bishop (1952) is chosen as a basis for further work. A further cor- 
rection to the Bishop approximation is developed to conform closely 
to experimental data. *- 

In Chapter III, an energy equation consistent with the assump- 
tions used in the elastic wave equation, coupled with the elastic 
equation of motion, is derived and is solved for the thermoelastic 
phase velocity and attenuation coefficient. 

Several interesting cases are solved numerically. These solu- 
tions indicate that the effects of thermoelastic damping are very 
important at extremely high frequencies. 



CHAPTER I 


WAVE. PROPAGATION IN AN INFINITE ELASTIC SOLID 
A. The Equations of Motion 

The equations of motion of an isotropic elastic medium are 
derived from basic principles in terms of cylindrical coordinates 
which were chosen in view of later applications to cylindrical rods 
Starting with the stress component equations of motion and with the 
strain components in terms of particle displacements, the displace" 
ment equations of motion are derived using Hooke's law. These equa 
tions are then shown to correspond to two types of elastic wave 
propagation in an unbounded solid. The first type consists of a 
volumetric expansion-contraction wave front called a dilatation 
wave. The particle motion in this wave is along the direction of 
propagation. The second type of wave is commonly referred to as a 
distortional wave. Particle motion in a plane distortional wave 
takes place perpendicular to the direction of propagation. In an 
unbounded solid, only these two types of wave are propagated. 

Figure 1 defines the cylindrical coordinate system along with 
the usual Cartesian coordinate system. The material of the rod is 
considered homogeneous and the section is in dynamic equilibrium 
under the influence of stresses and resulting accelerations. Dis- 
placements from the equilibrium condition are U along the radial 
coordinate r, V along the angular coordinate 0, and W along the 




FIG. 1. COORDINATE SYSTEM DEFINITION 



longitudinal coordinate z. The normal stresses are denoted by a ’ , 
ct q> and o z> and the shearing stresses are represented by T r g, t^, 
and t . The mass per unit volume is p. The stress equations of 
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motion for an elemental volume are 
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Adding the normal strain components yields the dilatation 
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Also useful in later derivations are the rotations to , ul, to 
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Within an isotropic, homogeneous material, each of the six com- 
ponents of stress at any point is a linear function of the six 
components of strain; 
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where A and (a are Lame's constants defined below. Poisson's ratio 
is denoted by v, and Young's modulus is E. 

* - laTrf < 6 > 

x “ a + via - 2v> ’ (7) 

The displacement equations of motion are derived from the stress 
equations of motion (1) by introducing the s tress -strain relations 
(5). After the appropriate substitution of various derivatives of 
the stress -strain equations and algebraic manipulations, an inter- 
mediate form of the equations of motion can be written: 
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The first terms of these equations can be recognized as the partial 
derivatives of the dilatation A with respect to r, 0, and z. The 
second terms can be shown to be the Laplacian of U, V, and W in 
cylindrical coordinates. The Laplacian in cylindrical coordinates is 


V s = 



+ - *r- 
r or 






Thus, the equations of motion (8) can be simplified: 


( 9 ) 


(A + n) ^ (u + 2 = pU (10a) 

L A * ^ |f + H^V + 4? (2 || - V) = pV (10b) 
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As a basis for further development, these equations are recast by 
using the derivation of the rigid body rotations u,., oig, and co z of 
equation (4) and by manipulating equation (8). The equations of 
motion in this form are known as the Pochammer-Chree equations: 
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B. Wave Propagation 


The Pochammer-Chree equations allow the determination of the 
propagation of two types of wave motion. 

Manipulation of equations (11) results in the wave equation 


a g A 

p at 2 " 


(A + 2|i) V 2 ^ 


( 12 ) 


in its most general form. This equation shows that the dilatation 
A is propagated throughout the elastic medium with a velocity 


Ci 


(A + 2u) 

P 



(1.3) 


More strictly, this wave should be termed irrotational. 

By differentiating (lib) with respect to z, and by subtracting 
a multiple of the differential of (11c) with respect to 0, the 
following form appears 

d 2 w 

M v 2 w . (14) 

ot^ p r 


Similar equations can likewise be obtained for Wq and Thus, it 

has been shown that the rotational wave form is propagated with a 
velocity 

C 2 =[n/p] l/2 . (15) 
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The strict term for this wave is equivoluminal; however, the 
more popular term is distortional or shear wave. 

It has thus been shown that, in an elastic solid, waves may be 
propagated with a velocity C x for waves involving no rotation and with 
velocity C 2 for waves involving no dilatation travel. The velocity 
designations, C-l and C 2 , are used for consistency with published 
literature . 

These results will be applied when the effects of the stress 
conditions at the solid boundary are introduced and the concept of 
dispersion is revealed. 
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CHAPTER II 


LONGITUDINAL ELASTIC WAVE PROPAGATION IN CYLINDRICAL RODS 

A. Elementary Theory 

It is possible, at least theoretically, to derive the propagation 
of stress waves in any bounded isotropic solid by solving the dis- 
placement equations of motion for the appropriate boundary conditions. 
In' fact, because of the reflection of plane elastic waves at free sur- 
faces and because of the extremely complicated fashion in which the 
dilatational and distortional waves are reflected as they travel up 
and down the bar, no exact solutions have been obtained. Thus, 
introducing boundary conditions into the equations of motion produces 
dispersion; that is, the velocity of wave propagation of a disturbance 
up and down the bar is a function of wavelength. 

The governing equations have been approximated by many in the 
field of elasticity. These approximations reduce to either of two 
types. The first introduces approximations into the equation of 
motion, while the second uses the exact equations of motion and only 
approximately satisfies the boundary conditions. The first method 
has been chosen for review and extension. 

Basically, three types of wave motion occur in thin rods; longi- 
tudinal, torsional, and lateral. In longitudinal motion, the elements 
of the rod extend and contract, but there is no lateral or transverse 
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motion of the axis of the rod. In torsional motion, each transverse 
section of the bar remains in its own plane and rotates about its 
center; again the axis of the rod is undisturbed. Lateral waves 
correspond to bending of the rod, and elements of the bar axis move 
laterally during the motion. 

Only longitudinal wave motion will be considered in this report. 
Historically speaking, the first treatment of longitudinal wave 
motion neglected the radial displacements arising from Poisson's 
effect. As a result, this first treatment predicted a wave velocity 
independent of frequency. A second elementary treatment was intro- 
duced by Lord Rayleigh, who shows that, by including the effect of 
radial inertia, the wave propagation velocity is a function of wave- 
length. These and further investigations, including the effects of 
radial shear by Mindlin and Herrmann and by Bishop, are compared in 
this chapter. 

The elementary theory of longitudinal wave propagation in a 
cylindrical rod assumes that the cross-sectional element faces remain 
plane and that the stress distribution is constant across the faces. 
Thus, the equation of motion can be directly written. Using Newton's 
law and introducing the material density p and forces due to the 
stresses on both sides of the element, we can write 


b 2 ^ _ ^ G z 


( 16 ) 
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By ignoring the lateral contraction due to longitudinal strain, we 
obtain 



Thus, the equation of motion can be written 


m E 

3t 2 p 8z 2 


(17) 


(18) 


Equation (18) is the classical first approximation for longi- 
tudinal wave motion in a rod. The velocity of propagation is 

C Q = [E/p] 1 / 2 . (19) 

The designation "C 0 " is used in the literature to distinguish it 
from the dilatational and shear velocities, C ! and C 2 . Notice, 
however, that the elementary treatment assumes that the plane trans- 
verse sections remain plane during the stress wave passage and that 
stress is uniform across the plane section. However, we know that 
longitudinal contractions and expansions will result in lateral 
expansions and contractions. The resulting lateral motion will 
cause non-uniform distribution of stress across the section of the 
bar, and plane transverse sections become distorted. The lateral 
effects will be small, however, when the wavelength of the wave is 
very much larger than the diameter of the bar. Thus we would expect 
any detailed analysis of longitudinal vibrations to reduce to the 
wave velocity C Q for very small ratios of the bar diameter to the 
wavelength. 
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B. Raylelgh-Love Correction 


The next important development in wave propagation in elastic 
rods was made by considering lateral contraction due to longitudinal 
strains. This was first published by Lord Rayleigh (John William 
Strutt) in 1877 and later by A. E. H. Love. Rayleigh's derivation 
was based on the principle that a frequency of free vibration of a 
conservative system can be obtained by equating the time average of 
the potential and kinetic energies of the motion. Love used 
Hamilton's principle in obtaining the equations of motion. The 
derivation used in this report will be developed from the displace- 
ment equations of motion. 

Because the rod is axisymmetric, the normal strain equations 
(2a), (2b), and (2c) simplify to: 



(20a) 


£ 


0 



(20b) 


£ 

Z 


3w 


(20c) 


The shear strain equations are all zero by virtue of axisymmetry and 
an assumed independence of z upon U. 

The radial strain e r simply is Poisson's ratio v times the 
longitudinal strain e z . Substituting equations (20a), (20b) and 
(20c) into the strain relations, and integrating with respect to r 
yield: 
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( 21 ) 


U = - vr M 
Bz 

Thus, equation (20b) becomes 



( 22 ) 


Now substitution of the above statements into the stress strain equa- 
tions (5) yields 


M - 2|xv g 

(23a) 

M - ,2nvg 

(23b) 

AA - 2p, g . 

(23c) 


These equations are then substituted into the radial equation of 
motion (la) which is then solved for the partial of radial stress with 
respect to r: 


Br 


vrp BzBt 5 


(24) 


Integration (24) with respect to r from 0 to r yields the radial 
stress: 


~ _ B^ r 2 
a r " vp BzBt 5 2 * 


(25) 
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Using equation (5) and the foregoing statements yields; 


E 



a 
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v 2 r 2 


dftt 

9 W ‘ 


( 26 ) 


Equation (23) is solved for the partial of longitudinal stress with 
respect to z, and is then integrated over the cross-sectional area of 
the bar to yield 


da 

z 



+ 


v 2 a g p d*W 
2 Sit 5 


Equation (27), substituted into the axial equation of motion 
results in the wave equation 


(27) 

(lc), 


E , v 2 a 2 d^ 

~pte s + ~ * ^5 • 


A solution of this partial differential equation is 


(28) 


w = e i.T}(.z-Ct) ^ 


(29) 


where T) - •£- , and L is the wavelength. The velocity of wave propaga- 
tion is then 


C = C, 


1 + 


p p p 

n v a 


^1/2 


(30) 


Thus, we see that the velocity of wave propagation is a function of 
the wavelength, known as "configurational dispersion." For long waves 
the wave speed is C Q , corresponding to the elementary theory. However, 


16 



equation (30) predicts zero wave velocity at zero wavelength. This 
is not borne out by experimental investigations, which show the 
wave velocities for short wavelengths to approach the surface wave 
velocity first discovered by Rayleigh. 

C. Mindlin and Herrmann's Method 

Mind 1 in and Herrmanns method for longitudinal wave propagation 
in a cylindrical rod proceeds from the displacements: 


U = “ U(z,t) 

(31a) 

V = 0 

(31b) 

W = W(z, t) , 

(31c) 


which are then substituted into an incremental energy equation 

= J dt J (T - V) dv, (32) 

t v 

and integrated over the cross section. The four resulting bar 
stresses are then removed by integrating the stress equations of 
motion over the radius of the section and by using Hooke's law to 
convert the stress equations of motion into displacement equations. 
To make the wave velocities fit the exact theory, the constants K 
and K x are then introduced into the displacement equations. The 
resulting displacement equations consist of two coupled partial dif- 
ferential equations: 
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a 2 K V - 8K^(A + n) U - 4aK 2 A ^ + 4aR = pa 2 


(33a) 


and 
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where 
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Mindlin and Herrmann assume the solutions to these equations to be: 


u = Ae i, ' (z ' ct) 
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As a result the dispersion equation was written as 
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A + 2u P 

8K?A 2 

P a n p J 

” L. 
P 

“ 2 2 2 
a 2 r) 2 p 2 


= 0 . 


(35) 


This equation represents two modes of wave propagation in the bar. 
For long wavelengths, the wave velocities are 



(36a) 


and 


C 2 00. 


(36b) 



The first branch (equation (36a)) agrees with the elementary theory. 
For short wavelengths, the wave velocities are 


C 2 ■ K a Ut 
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(37a) 


and 


,2 _ A + 2u 


(37b) 


The first parameter K 2 is used to fit the dispersion curve to the 
Rayleigh surface wave velocity. The second parameter K 2 is used to 
fit the dispersion curve to the Pochammer-Chree curve at a common 
point C 2 = 2c| for all values of Poisson's ratio v. The coupled 
equations of motion (33a) and (33b) are then written to eliminate U. 
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(38) 


Inspection of equation (38) shows that Mind 1 in and Herrmann's method 
improves the Rayleigh-Love differential equation by adding the radial 
shear effect and by improving the radial inertia term. 

D. Bishop's Method 


Bishop's original paper on the velocity of wave propagation in 
circular elastic rods was prepared in 1952 in an effort to improve 
the prediction of elastic wave propagation velocity as a function 
of wavelength. Bishop derived an approximation to the equations of 
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motion in Cartesian coordinates including the lateral inertia as 
Rayleigh and Love and also an approximation to the lateral shear 
stress. He used a theorem relating to convex boundaries to remove 
geometric problems in describing the cross sectional boundary. How- 
ever, for clarity, we shall use Hamilton's principle and cylindrical 
coordinates to derive the equation of motion. 

The strain relations for e r , €q, and e g are the same as that of 
Rayleigh and Love. Unlike Rayleigh and Love, Bishop did not ignore 
the shear strain. 
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rz 


Su 

dz 
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r IbP • 


(39) 


Hamilton's principle states that the first variation of the time 
integral of kinetic energy minus the potential energy during a process 
is zero, as symbolically stated below: 



(T - V) dt = 0. 


(40) 


The strain potential energy of an elastic material (per unit volume) 
is the integral of the elastic force over the displacement 


V = 


F dz 


(41) 


where 


F = Ez. 


(42) 
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Integrating from zero to e z results in the strain energy per 
volume of 

V = 1/2 Eef. 

Likewise, the shearing strain energy is 

V = 1/2 :W £, 

Thus, the total potential energy is 


r 9 z 

v = 1/2 J* J J (Ee 2 + h 7 2 z ) rdrdedz, 
0 0 0 
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or, for the rod in question, 


r z 

V - it J J' (Ee 2 + M7r Z ) rdrdz. 
0 0 


Substitution of the strain relations yields 
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The kinetic energy per unit volume is 
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Substituting the strain relations and integrating, we obtain 
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0 0 


When we substitute these forms into Hamilton's equation and integrate 
over the cross section of the bar, we obtain (introducing subscript 
notation for partial derivatives) 


tg 2, 
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ita‘ 
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Forming the functional 
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and using the Euler-Lagrange necessary condition of the variational 
calculus, we obtain 


H - g Of/3w 2 ) - g (SF/aw t ) + g* (»/aw zz ) + gs (SF/3w tt ) 
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zt 


(52) 


After we form and substitute the indicated partial derivatives, the 
resulting equation of motion is 
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By noting the elementary propagation velocity, C Q , and the distor- 
tional wave velocity in an extended medium, C 2 , we can simplify 
equation (53): 
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Again, assuming a solution 

w = 

yields the dispersion equation 


(55) 
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(56) 


For long wavelengths (t] approaches zero), the wave velocity is the 
same as that for the elementary theory. For short wavelengths (r) 
approaches infinity) , the velocity of wave propagation approaches 
the velocity of shear (dis tortional) waves in an infinite medium. 
Thus, we observe, by comparing equation (56) with equation (30), 
that Bishop's contribution is the addition of the radial shear term 
to the Rayleigh-Love analysis. 
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E. Comparison of Bishop's Method and Mindlin and Herrmann’s Method 




Comparison of the results of these two methods begins by compar- 
ing their dispersion curves with numerical approximations to the 
Pochammer-Chree equations by Bancroft and Hudson in figure 2. It 
can be noted that both the Bishop and the Herrmann and Mindlin dis- 
persion equations result in wave velocities greater than the 
Pochammer-Chree numerical solutions for all wavelengths. However, 
the errors resulting from the Mindlin and Herrmann method are generally 
smaller, due in fact to the use of the constants K and K x for curve 
fitting. For our purpose an insight into the complexity of the 
Mindlin and Herrmann method can be gained from a comparison of the 
equations of motion of the methods. 

Using Mindlin and Herrmann’s one-dimensional equation of motion 
(38) and assuming that 


W = 

(57a) 

and 


R = Z = 0, 

(57b) 

we obtain the following form: 


- W " W + (A + 2(j.) + K p, W 

P ZZ tt 8 r 2 (a + ^) L J 


pa 2 w a% 2 u(A + 2u) w _ 0 

8K 2 (A + jo.) ttt:t 8K 2 (A+|i)p ZZZZ 

(58) 
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which is to be compared with Bishop's equation of motion (53) 


E tt 4- vg a 2 TT 

P W 2Z • "tt + “I - W ttZ2 


2p zzzz 


0. 


(59) 


F. Improved Longitudinal Wave Equation 

The simple form of the Bishop equation of motion suggests its 
use in an improved form as an engineering tool. The basic equations 
of motion are reformulated here in terms of the stress equations with 
constants K and Kj. used to fit the resulting dispersion curve to the 
numerical solution of the Pochammer-Chree equations. 

Starting with the stress-strain relationships of the Rayleigh- 
Love analysis (equation (23)) and the shear strain relation of the 
Bishop analysis (equation (39)) 


a r - ^ - 2,v | 

(60a) 

a 0 = M - 2pv g 

(60b) 

a 2 = M - 2pv ^ 

(60c) 

- - 

T - - |J.vr -r — 5 

rz n dz 

(60d) 


and substituting into the radial stress equations of motion (equa- 
tion (la)), we obtain the following: 


^2 

Sr 


= K 2 |_ivr 


S 3 W 

Sz 3 


- K‘ 


S 3 W 

vrp SzSP 


(61) 


V 
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Integrating equation (61) with respect to r yields 


r _ K 2 uvr 2 d 3 W K x v P r2 ^ 

a r 2 ^ + 2 ~~ iF * 


( 62 ) 


Now using equation (59c) and the strain equation (2) and taking the 
derivative with respect to z yields 



After we integrate over the cross section and divide by the cross- 
sectional area to find the total bar stress, and when we substitute 
into the classical bar stress equation of motion (equation (16)), 
we have the improved wave equation written with previous notations; 



k 2 c| 


v 2 a 2 b\ . „2 v 2 a 2 bSl b 2 ^ 

— bP + K 1 — • aF = 0 


(64) 


Introducing the solution 


W = e iT l (z " Ct) 


(65) 


yields the dispersion equation 



+ C 2 


1 + 


iq v 2 a 2 if 


2 2 P 

2p 2 v a 2 n _ 
2 


K 2 C 


0. 


( 66 ) 


For long wavelengths (t) ->0), the wave speed equals the elementary 
wave speed, C 0 . For short wavelengths, the wave propagation velocity 
is 
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( 67 ) 



The ratio of the constants can be used to fit the dispersion 
curve to the Rayleigh surface wave velocity. This ratio equals 92.6 
percent of the shear (distortional wave velocity in an infinite 
medium for a Poisson's ratio of 0.29, which is typical of steel. 
Then, if the dispersion curve is fit at a/L ■ 0.293 and C/C 2 ■ 2, as 
in Herrmann and Mindlin, the results can be plotted (see figure 2). 
Also shown for comparison are the Bishop curve and Mindlin and 
Herrmann curve (Pochammer- Chree numerical solution by Hudson). 
Notice that the improved solution lies between the exact and Mindlin 
and Herrmann curves. Thus, we have a very simple improvement to the 
Bishop form of the longitudinal wave equation. 
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CHAPTER III 


THERMOELASTIC WAVES 
A. Thermoelastic Equations of Motion 

Chapter I reviews the basic principles of elastic wave propaga- 
tion in an unbounded solid and Chapter II examines the characteristics 
of longitudinal wave propagation in cylindrical bars. This chapter 
discusses the salient points of thermal modification of longitudinal 
elastic waves. 

When a longitudinal wave passes through a bar, any given element 
is compressed and extended. These volumetric changes are accompanied 
by heating and cooling. The heat generated during the compressional 
phase is conducted through the bar. The resulting temperature change 
affects the state of strain by the coefficient of thermal expansion. 

The approach taken here is to combine these thermal effects with 
the improved wave equation developed in Chapter II. The resultant 
effects upon the purely elastic wave motion are discussed in terms of 
a change in the propagation velocity, and the addition of a damping 
effect upon the mechanical energy of; the elastic wave due to heat 
conduction. 

The thermally coupled case of longitudinal harmonic wave propa- 
gation is derived from the stress-displacement radial equation of 
motion developed for the improved elastic wave equation ( 62 ) coupled 
with the energy equation. 
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The thermoelastic stress-strain relationship is written: 


€ r - 5 « T +| [a r - v(cr e + <I Z )] (68) 

where Oi - coefficient of volumetric thermal expansion. By observing 
that a r equals Oq, we write 


Scr 

2 

Sz 


= E 


S 2 ^ 


§* " K2 ^ v2r2 S* + spS* ' T s 


bSl 


2 p p 

1 1^- r-.V» ^ 


S^W 


eq: St 


(69) 


This equation is then integrated over the circular cross section to 
yield: 



E 


S 2 W 

Sz 2 


K 2 


a 2 uv 2 S*W . ..a a 2 v 2 o 2 S'H? Ea St 
2 + K i 2 Sz 2 St 2 3 Sz * 


(70) 


Equation (68) is introduced into the stress equation of motion derived 
for the elementary case, 


*k. 

Sz ^ St 2 ’ 


(71) 


and, thus, the following equation of motion is obtained: 


— w - 

p zz 2 


K x v 2 a 2 Ea _ 

+ 0 W . . - o t — w ... 

zzzz 2 zztt Jp z tt 


(72) 


The general form of the energy equation is written most often in the 
form 


ggo S(div U) 
3(1 - 2v) St 


+ pC 


St 
e St 


= K^T + Q, 


(73) 
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where 


C £ * specific heat at constant strain 
U * displacement vector 
Q = internal heat addition 

and 

k = thermal conductivity. 

But div U is the dilatation A, which is evaluated to be 

A = (1 - 2v) H . (74) 

The Laplacian of temperature, is evaluated by assuming that the 

temperature gradient across the cross section is zero. Thus, the 
energy equation (with no internal heat generation) becomes 


St 


, E aT rT 
+ — r — W . 
3 zt 


k 

P 


T . 
zz 


(75) 


Equations (72) and (75) are the coupled partial differential equa- 
tions of thermoelastic disturbances in a cylindrical bar of infinite 
length. Because we wish to find thermoelastic damping attenuation 
coefficients, as well as phase velocity, solutions of the type 


and 


W = W e 
o 


i(TjZ-.Wt) 


T = T e i( ^ Z “ wt) 
o 


(76a) 

(76b) 


are most fitting. In general u and t) are complex quantities. The 
wavelength of the plane harmonic waveform is 2jt/uteal* 
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Substituting equations (76) into (72) and (75) yields the two 


relations: 


W„ 


E „ K 2 v 2 a 2 u . , K l v2a2 
p ’> 2p""* ’l + ~2— 


*1 2 U 2 + U 2 


O' 


EOT 

c 

3p 


ir) 


(77) 


and 


>T 0 W 0 Ea 
— ^^o 




(78) 


Eliminating T and W and multiplying by -i/C e yields the thermoelastic 
wave equation for longitudinal harmonic plane waves propagating in an 
infinite cylindrical rod 



If w is chosen as a real quantity representing waves of assigned 
frequency, the complex wave number q can be found as a function of 
frequency. By expanding the form in equation (76), we get 
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e 


lma 8 2 cos (z - -J2— t) . 

' \eal / 


e i(T)z-wt) 


1 [ Val + 1 ’W'“ t] ^ 

“ • 6 


( 82 ) 

Thus, the wave velocity equals w/ilreaij and the attenuation coefficient 
q equals T)^ tna g. Thus, the damping is a function of wave travel and is 
independent of time. Solutions of assigned wavelength can also be 
obtained from equation (79) but were not found due to the general 
interest of the assigned frequency solutions. 

B. Uncoupled Solution 


The thermoelastic wave equation above consists of a pure thermal 
wave coupled through the constant e to the improved elastic wave equa- 
tion found in the last chapter. Thus, when the coupling constant is 
zero, the thermal and elastic waves exist separately. If equation 
(79) is uncoupled (e = 0), the solutions are 


ipC 


'imag 


w, 


(83) 


corresponding to the thermal mode, and 
K 2 v 2 a 2 w 2 \ 


•real 


l-i./ 2 


C 2 - -i 

L o 2 


K 2 v % 2 




-1 + 1 + 


2K 2 v 2 a 2 (j,w 2 


i 1 / 2 


2 Kfv 2 a £ 

p l c o *T 




(84) 
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corresponding to the elastic mode. The phase velocity, V « W / T l r eal’ 
is solved in terms of the wavelength L ■ 2rc/w and yields 


1 + t 


3 


1 + 


2K 1 2 v 2 a 2 jt 2 


v a « 
L 2 _ 


(85) 


This equation is identical to the dispersion curve relationship for 
the improved elastic wave in Chapter II and represents an imaginary 
material for which a - 0. However, a ■ 0 is not a valid physical 
approximation as we shall see in the following discussion. 


C. Coupled Solution 


For all materials considered, the coupling coefficient is very 
small. This fact, along with the fact that a computerized solution 
of the sixth order roots of equation (79) would be required, sug- 
gested the use of a perturbation technique. To condense the notation, 
we use a characteristic frequency co* to define the nondimens ional 
quantities as suggested by Chadwick (1962): 


0J 


1 w* 

pC C 2 
cJ* 


( 86 ) 

(87) 

( 88 ) 
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Using these quantities defined in equation (79), we obtain 


(X + i| 2 )[| 2 (l - d 2 v 2 X 2 ) + c 2 v 2 g 4 - X 2 ] + e£ 2 X - 0, (89) 


where 


K 2 a 2 |ap 2 C 2 C 2 

° 2 “ 2 » 
2Ek 2 


and 


(90) 


d 2 


2 ” 2 p 2 C 2 C 2 
H e o 


K ? 


Ik 2 ” 


(91) 


When X is regarded as a real constant (assigned frequency), equation 
(89) is of the fourth order in g. However, in this discussion, the 
particular root corresponding to the quasi-elastic wave is desired, 
and therefore the form of £ chosen is 


t = Ip + 5g, (92) 

where £ represents the improved elastic wave solution when e = 0. 
The solution of equation (88) under this condition is 


-(1 - d 2 v 2 X 2 ) + 4 (1 - d 2 v 2 X 2 ) + 4c 2 v 2 X 2 
2c 2 v 2 


(93) 


The positive sign was chosen by comparison with equations (84) and 
(85). 


When equation (92) is substituted into equation (89) and if 
higher order terms are neglected, the result can be expressed as 
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* 


(94a) 


-2e| 0 X a O l/2 + e) + 12€|^ 1 / 2 

£ * | Q 4* " 

4 X 2 (p 1 / 2 + €) 2 + 4|^p 2 

where 

(1 - d 2 v 2 X 2 ) 2 + 4c 2 v 2 X 2 = p. (94b) 

Equation (94) is of the form (when compared to equation (80) 


I = 


c o<! + l 30- 


CO" 


(95) 


Therefore, the quasi-elastic velocity V is 

C X 4C 0 X[X 2 (p l/2 + e) 2 + l*P 2 ] 

V » T- 2 — = 2 > 

5 real -g £l - 2eX 2 (p 1 / 2 + e) ] 

and the damping coefficient q is 


(96) 



eco *! 3 p 1 / 2 
b o K 

2C 0 [X 2 (p l/2 + e) + l^p 2 ] 


(97) 


Equations (96) and (97) represent the main contribution of this 
report to wave motion in a bounded thermoelastic solid. 


D. . Numerical Results 

When equations (91) and (94) are solved numerically for the 
ratio of the thermoelastic velocity to pure elastic velocity (t) 
as a function of the reduced frequency X, the results may be plotted 
as shown in figure 3. The data throughout this section are typical 
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of structural aluminum and are shown for bar radii between 10” 3 and 
10 2 centimeters. Figure 3 makes clear that thermal coupling modifies 
the elastic wave velocity by a small amount at low values of X. At 
higher values of X, the modification becomes negligible. The fre- 
quency at which the thermoelastic wave velocity asymptotically 
approaches the pure elastic velocity is a function of the bar radius. 
The characteristics of these curves are the same as in Chadwick. The 
relationship of bar radius can be eliminated from figure 3 by changing 
the independent variable to the ratio of bar radius to wavelength 
(a/L) (see figure 4). Thus, figure 4 shows that the thermoelastic 
modification of elastic wave propagation velocity is primarily depend- 
ent upon the wavelength in relation to the bar radius. 

Figure 5 presents the behavior of the thermoelastic damping coef- 
ficient (8) as a function of the reduced frequency (X) with bar radius 
as a parameter. The damping coefficient (q) varies as the square of 
the reduced frequency at low frequencies and asymptotically approaches 
a finite limiting value that is inversely dependent upon the square 
of the bar radius. It is clear then that very short waves traveling 
in small diameter rods are severely attenuated. When the same data 
are plotted as a function of the ratio of radius to wavelength (fig- 
ure 6), the breakover point of the curves occurs when the wavelength 
is roughly equal to one half of the bar diameter. 

The following is an example of the magnitude of this damping: 
a 10" 3 cm wave traveling in a bar of the same radius would be damped 
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to half amplitude in 6.93 cm. As a more practical case, a one centi- 
meter wave in a bar of the same radius would decay to half amplitude 
in 69.3 kilometers. 

E. Limiting Case for a/L Approaching Zero 

The above results for thermoelastic rods can be compared with 
Chadwick's results for an infinite thermoelastic solid by noting 
the behavior of equations (93), (96), and (97) at large wavelengths. 
When a/L approaches zero, equations (90) and (91) approach zero and 
equation (89) becomes 


(X + i|) (| 2 - X 2 ) + e| 2 X = 0, (98) 


which is exactly Chadwick's result. 

The uncoupled solution of equation (98) is 


I 


2 

O 


= X 2 


(99) 


and 

tl - iX. (100) 

Equation (99) corresponds to the elastic wave and equation (100) to 
the thermal mode. Using equation (99) and the perturbation form, 
equation (92), and ignoring higher orders of e, we obtain 


and 


6 = X(1 - O 

Wl 2[(1 + 2e) + X a ] 


( 101 ) 


I 


eX 2 

imag 2 [ (1 + 2e) + X 2 ] * 


( 102 ) 
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When equation (101) is used with (95), the thermoelastic velocity- 
can be written as 


V '°0 1+ wfx?)- (103) 

This form for the wave velocity can be shown to correspond with 
Chadwick's results of the first order in e. At large wavelengths 
(X — > 0) , equation (103) reduces to 

v x^0 ■ V 1 +!>• < 104 > 

This result is easily seen in figures 3 and 4 if we realize that e 
for aluminum is 0.0054. Likewise, equations (95) and (103) are used 
to solve for the damping coefficient, 


q 


eX* 


2(1 + X 2 ) 


(105) 


Equation (105) corresponds to the results of Chadwick (to the first 
order in e) and is shown in figure 5. The reason that the damping 
coefficients predicted by the present theory exceed the value 
obtained for the limiting case may be seen by taking the ratio of 
equation (97) to equation (105). At small wavelength (X -»<»), equa- 
tion (105) reduces to 


'X 00 


0j' ? £ 

2C 


(106) 
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The elastic velocity of the present solution is C Q instead of C x in 
the Chadwick solution because, in the present theory, the dilatation 
is expressed as 

A “ (1 " 2v) || (107) 

instead of the form of equation (3). The thermoelastic behavior 
for small values of radius to wavelength ratio is the same, with the 
exception of dilatational velocity, as reported in Chadwick's analysis 
of an infinite thermoelastic solid. 


F. Thermoelastic Bar Number 

The general behavior of the damping coefficient can be more 
clearly seen if we observe that, in practice, X will be very much 
smaller than one, since i o* is in the order of 10 11 radians per 
second. Using this fact, we can approximate the behavior of q as a 
function of X, Consider the behavior of equations (93) and (94) when 
dvX (thus, also cvX) is less than one. Equation (93) yields 


| = X, 

b o ' 


and equation (97) yields 


q = 


to*eX 2 

2C o * 


(108) 


(109) 
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Equation (105) neatly expresses the behavior of the damping coeffi- 
cient as a function of X 2 . When dvX is greater than one, equation 
(93) yields 


*o " c X • 


( 110 ) 


and equation (97) yields 


_ t/ed 

q 2C 0 c 3 v S * 


( 111 ) 


Equation (107) then predicts the value of q at very high frequencies. 
At the particular value of X for which dvX is identically one, equa- 
tion (93) reduces to 


i -x. 

b o cv 


( 112 ) 


The damping coefficient, equation (97), is thus 


q = 


w*e 


4c C 2 v 2 "7 
o d 


(113) 


This point corresponds to a maximum slope condition and can be used 
as a demarcation point for the approximations contained in equations 
(109) and (111). Thus, the behavior of the damping coefficient can 
be approximated by equations (109), (111), and (113). A specific 
example is shown in figure 7 for an aluminum bar of 1 centimeter 
radius. Thus, the parameter dvX has been shown to be of primary 
importance to the behavior of the thermoelastic damping coefficient 
in a cylindrical rod. 
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CONCLUSIONS 


This study considers longitudinal thermoelastic wave propagation 
in cylindrical rods. More specifically, the thermoelastic effects on 
an assigned frequency wave are shown to be small in terms of the 
propagation velocity. However, in comparison, the thermoelastic 
damping effect is large for very high frequency waves traveling in 
small diameter bars. The behavior of the thermoelastic damping coef- 
ficient is linked to the "thermoelastic bar number" dvX, developed 
herein. Finally, approximations are developed, based upon the value 
of the thermoelastic bar number, to predict the behavior of the 
thermoelastic damping coefficient as a function of frequency. 

It is hoped that this 'report provides a basis upon which 
investigations into thermoelastic effects on other structural 
elements will proceed. 
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